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Abstract. - We investigate quantum effects in the mechanical properties of elastic beams on 
the nanoscale. Transverse quantum and thermal fluctuations and the nonlinear excitation en- 
ergies are calculated for beams compressed in longitudinal direction. Near the Euler instability, 
the system is described by a one dimensional Ginzburg-Landau model where the order param- 
eter is the amplitude of the buckling mode. We show that in single wall carbon nanotubes with 
lengths of order or smaller than 100 nm zero point fluctuations are accessible and discuss the 
possibility of observing macroscopic quantum coherence in nanobeams near the critical strain. 



Introduction. - The progress in miniaturization of electromechanical devices towards the 
nanometer scale (NEMS) is beginning to reach the limit, where quantum effects play an im- 
portant role PJElEl- For example, in nanoscale beams phonons may propagate ballistically, 
leading to a quantized thermal conductance [Ij . Moreover a sizeable contribution to the forces 
between plates and beams which are separated by less than one micron is the Casimir force 
between neutral objects due to the modification of the electromagnetic vacuum The 
combination of electrical and mechanical properties may be studied via quantized transverse 
deflection due to charge quantization of charged, suspended beams in an electric field (7j. Sim- 
ilarly the standard Coulomb-blockade in small metallic islands or in semiconducting quantum 
dots may be used to mechanically transfer single electrons with a nanomechanical oscilla- 
tor jSJE]. Regarding possible applications of nanomechanical sensors, Si-based resonators in 
the radio- frequency regime were recently fabricated and manipulated . In the present work 
we focus on quantum effects in mechanical resonators on the nanometer scale, in particular 
in single wall carbon nanotubes (SWNT). Due to their small masses and remarkable elastic 
properties down to nanometer scale, carbon nanotubes are ideally suited to study effects like 
phonon quantization 11, the generation of non-classical states of mechanical motion |12) 
or macroscopic quantum tunnelling out of a metastable configuration J3j. On the classical 
level, both the thermal Brownian motion of single nanotubes clamped on one side jjll and the 
discrete eigenmodes of charged multiwall nanotubes excited by an ac-voltage jj^l have been 
detected experimentally. More recently, the thermal vibrations of doubly clamped SWNT's 
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down to lengths of around 0.5/im have been observed with a scanning electron microscope 
In all of these cases it turns out that the measured transverse vibrations of nanotubes agree 
reasonably well with the predictions of an elastic continuum model. Its applicability even on 
the nm scale is also supported by molecular dynamics simulations which show that SWNT's 
down to lengths of around lOnm are well described by an effective elastic continuum, respond- 
ing in a reversible manner up to large deformations |17j . In the following, we will therefore 
use the standard theory of an elastic continuum | I18| for carbon nanotubes which are clamped 
between two fixed end points. We calculate both thermal and quantum fluctuations of the 
nanotube under longitudinal compression, including properly the nonlinearity in the bending 
energy. It is shown that in SWNT's with a length below lOOnm the crossover from thermal 
to quantum zero point fluctuations is reached at accessible temperatures of around 30mK. 
We also discuss the possibility to realize coherent superpositions of macroscopically distinct 
states by observing the avoided level crossing near the degenerate situation above the critical 
force of the well known Euler-buckling instability. 

The Model. - Our model system is a freely suspended SWNT of length L and diameter 
D which is fixed at both ends, allowing only transverse vibrations. In addition we consider a 
mechanical force F which acts on the beam in longitudinal direction (F > for compression). 
In a classical description the beam is then completely described by the transverse deflection 
(/)(s) parametrized by the arclength s S [0, L]. We assume the beam to be incompressible 
in longitudinal direction and only keep a single transverse degree of freedom for simplicity 
(see below). For arbitrary strong deflections 4'{s) the nonlinear Lagrangian of the system is 
then [El 



(1) 



Here a = m/L is the mass density, while the bending rigidity /i = EI is the product of the 
elasticity modulus E and the moment of inertia / = TrD^d/8, with d an effective wall thickness. 
For small deformations |0'(s)| <ti 1 the Lagrangian is quadratic, leading to the standard linear 
equation of motion 

a4> + + = (2) 

for the transverse vibrations of an elastic beam under compression. The corresponding eigen- 
modes (pn and eigenfrequencies Un depend on the boundary conditions. We assume that the 
experimental reahzation is well described by clamped ends at both sides, 0(0) = 4'{L) = 
and 0'(O) = (f>'{L) = 0. The exact (f>n's are then given by a superposition of trigonometric and 
hyperbolic functions, and the a;„'s by solving a transcendental equation , ,18) . In the follow- 
ing, for some of the analytic expressions, we will use boundary conditions without bending 
moments at the ends of the beam, 0"(O) = 4'"{L) = 0. This leaves the essential physics 
unchanged and permits one to write down simple expressions for the eigenfunctions in the 
normal mode expansion 

0(s) = sin — s (3) 



and its eigenfrequencies 



— ■ (4) 



Clearly the modes soften with increasing compression F, up to a critical force Fc = //(tt/L)^ 
where the fundamental frequency a;„=i(F) vanishes. Then the system reaches a bifurcation 
point, the well known Euler instability, beyond which 0i ~ sin (jrs/L) becomes the new stable 
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Fig. 1 - Mean square displacement at the center of a clamped beam in linearized theory for the free 
{S = 1) and compressed {S = 1/4) case. The dotted lines show the classical behaviour if one considers 
only the leading contribution of the first mode, the dashed-dotted line indicates the increasing zero 
point fluctuations versus decreasing crossover- temperature, {a'^{5,T 0),T*{5)}, for compression 
approaching the critical force. 

solution of the static problem. For clamped boundary conditions with finite bending moments 
at s = 0, L the critical force is four times larger, and the shape of the stable solution in the 
static problem for F > Fc has the form sin^ (ns/L). Near criticality, the frequencies of higher 
modes n = 2, 3, . . . remain finite. The dynamics at low frequencies is thus determined by the 
fundamental mode alone. The nonlinear field theory eq.(^ may be quantized in the standard 
manner by requiring canonical commutation relations [4>{s,t),JI{s' ,t)] = ihS{s — s') between 
the field (j) and its canonically conjugate momentum 11 — a<j) at equal times. In the linear 
regime, the problem is reduced to an infinite number of harmonic oscillators. Introducing 
oscillator lengths I'j. = h/{mkiOk) with k = mr/ L and n = 1, 2, . . . , the amplitudes 

Ak = ^{al + ak) (5) 

are expressed by the standard creation and annihilation operators aj^ and Uk ■ The effective 
masses arising in Ik turn out to be rrics ~ 3/8 of the beam mass for the fundamental mode 
but are generally mode dependent for clamped boundary conditions. 

Thermal vibrations. - In the linearized theory, the mean square displacement of the 
beam is trivially calculated from the normal mode expansion eq.( Assuming a thermal 
occupation of the discrete phonon modes one obtains a maximum value at the center of the 
beam, which for undamped boundary conditions reads 

n odd V \ / 

1 

Here the scale is set by the oscillator length Iq = (h/rricffUJo)^ and the temperature Tg = 
huJo/kB is associated with the frequency scale of the fundamental mode loq = uJi{F 
Fc)S~^ ~ 1.02uji{F — 0). The parameter S = {Fc — F)/Fc determines the dimensionless 
distance from the critical compression force. At temperatures larger than Tq one obtains the 
usual equipartition theorem result, where cr^ ^ T/ujq increases linearly with T as observed 
on sufficiently long nanotubes |14II16| . For low temperatures the mean square displacement 
remains finite due to zero point fiuctuations. As shown in Fig. ^ the crossover to this regime 
occurs at T* ~ OATq in the absence of an external force 6=1, giving accessible temperatures 
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of around 30mK for typical SWNT's (see Table QJ. Unfortunately, the associated transverse 
deflection amplitude Iq is only of order lO^^nm and thus beyond accessibility of standard 
displacement detection techniques. There are a number of ways, however, to measure such 
tiny amplitudes, for instance by capacitively coupling the beam to the gate of a single electron 
transistor |2()j or using its electrostatic interaction with a free-standing quantum dot |2H . Of 
course the fluctuations are strongly enhanced near the critical force Fc, where only the fun- 
damental mode n — 1 contributes and thus (T^(r — 0) increases like Iq/\/S. In this case, the 
problem of measuring the thermal to quantum crossover is somehow inverted, since now the 
decreasing crossover temperature T*{S) ^ TqVS poses a limiting factor for an observation. 
Moreover the divergence in eq.(|SJl for S marks the breakdown of the linearized theory. 
To explore the quantitative enhancement of the fluctuations and the change of the relevant 
length and energy scales near criticality, one has to include the nonlinear terms in the bending 
energy eq.(nj. 

Buckling instability. - In order to describe the behaviour near the buckling instability, we 
insert the Fourier expansion eq.(|2Il into the corresponding nonlinear Hamiltonian. Keeping 
only the fundamental mode, since all higher modes have no influence near criticality due to 
their nonvanishing frequencies, the interacting field theory is reduced to a one particle problem 
in terms of the coordinate Ai, and its canonically conjugate momentum V = —ihd/dAi. 
The force term generates a negative contribution to the quartic term in Ai, driving the 
system unstable. The nonlinearity in the curvature term, however, over-compensates this and 
guarantees stability even for fixed length of the beam jJS] . For clamped boundary conditions 
one can use the approximate shape sin^ (tts/L) which becomes exact near (5 — > 0. One ends up 
with a quantum mechanical one particle Hamiltonian with an anharmonic oscillator potential 

^ = 2^ + l^cffC^oM? + ^-fAt (7) 

and an anharmonic coefficient 64 — {tt/ L)^FcL. A similar effective description of quantum 
effects near the buckling instability has been derived in The nonlinearity there, however, 
arises from longitudinal stretching while we keep the length of the nanotube fixed. It is now 
convenient to define a dimensionlcss coordinate yhyAi = ly, where I — Iq (27r^)~^/^ (L/Iq)^^^ 
is the characteristic magnitude of the deflection, where the quartic term due to the nonlin- 
ear bending energy is of the same order than the kinetic energy. The Hamiltonian is thus 
transformed to a dimensionlcss form 

H = hLO f-V + i£y2 i 4 
V 2 2(5'^ 4^ 

with u) — h/ (nicffP) as the characteristic frequency scale near the critical compression force 
Fc- It differs from the fundamental frequency loo of the classical transverse vibrations by a 
factor 

which also determines the size 6 of the critical regime. For SWNTs with length L ~ 0.1/xm, 
S^^'^ is of order or smaller than 10~^ (see TableOJ. The potential energy in eq.(|Sl) exhibits the 
standard Landau bifurcation from a single to a double well as the external force is increased 
through its critical value Ff.. Indeed our zero dimensional quantum problem is equivalent 
to a one dimensional classical Ginzburg Landau theory |23| . Let us consider first the mean 
square displacement at the center of the beam. In the harmonic approximation this diverges 
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Fig. 2 - Meansquare displacement at zero temperature at the beams center near the Euler instabihty. 
The solid line is the numerically exact cr^ with mean zero, while the dashed line shows the harmonic 
approximation for F < Fc and the dotted line the result ~ J/min corresponding to the new stable 
minimum in the broken symmetry phase F > Fc- 

Fig. 3 - Numerically exact excitation frequency (solid line) and tunnelsplitting A (dashed line) ac- 
cording to ea.( null in units of oj near the buckling instability. The dotted line corresponds to the 
harmonic approximation for F < F^. 



as F approaches the critical value from below. For F much larger than Fc it is simply 

~ 1/2 

determined by the stable minimum of the effective Landau energy at j/min — ± (I'^l/'^) 
giving CT^ — P\5\/5. As shown in Fig. [3 the exact result smoothly interpolates between 
those two limits giving a finite value (t(Fc) = 0.68 Z at F^ which is of order O.lnm for typical 
SWNT's (see Table Pi . A similar behaviour is found for the lowest excitation frequency of 

~ 1/2 

the beam. In the harmonic approximation it vanishes like uJi{F) = w • ((5/(5) . Above the 
critical force the lowest excitation is the small oscillation in one of the degenerate minima of 
the anharmonic oscillator. This is true, however, only in a classical description. Quantum 
mechanically, the lowest excitation is the exponentially small tunnelsplitting A which lifts the 
degeneracy between the two states localized in the left or right well of the effective potential. 
Again, the exact numerical result for uj = {Ei — EQ)/h starts to deviate from the harmonic 
expression at around 6 5 6 and approaches a finite excitation frequency to — l.lu) aA, Fc 
(see Fig. Ol. For S < —3(5 it vanishes exponentially in good agreement with the WKB result 
ea.( llOf) for the tunnelsplitting A. It is remarkable that the excitation frequency precisely 
at Fc, which is of order 2tt ■ O.OlGHz for the parameters of Tabled is no longer related to 
the characteristic frequency ujo of the classical problem but scales like w = uJo{4:huJo / FcLy , 
remaining finite only through a genuine quantum effect. Unfortunately, the smallness of the 
size 6 « 10^* of the critical regime requires fine tuning the compression force F very close 
(S ~ (5) to its critical value in order to see deviations from classical behaviour near the buckling 
instability. 

Macroscopic Quantum Coherence (MQC). - In the regime beyond the critical buckling 
force, the state of lowest energy corresponds to a stationary finite deflection amplitude ymin = 
i (l^l/<^) which is lower in energy by ?iw ((5/(5)^/4 than the configuration with no deflection. 
The direction in which the buckling occurs is arbitrary, however, in a realistic setup like that 
in ref. ^BJ the boundary conditions in the buckled state are likely to break the perfect rotation 
symmetry assumed in |23]. In this situation, the transverse deflection is described by a single 
degree of freedom with only two degenerate states. Quantum mechanically, these states are 
split into a narrow doublet with energy separation fi/S. due to tunnelling. Sufficiently far above 
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Table I - Characteristic parameters for quantum effects in SWNT's of length L = O.l/^m and 
diameter D = 1.4nm. We assume a Young's modulus E — ITPa and an effective wall thickness 
d = 5 ■ 10~^nm. These parameters are consistent with recent measurements (16). 





JTloff [kg] 


wo/27r[GHz] 


lo [nm] 


ro[mK] 


Sl/2 


Ld/2n[MUz] 


Z[nm] 


0.19 


1.3 • 10"^^ 


1.4 


0.01 


65 


0.6 • 10"^ 


8.5 


0.13 



the critical force its magnitude may be determined from a WKB-calculation in the anharmonic 
oscillator potential, giving 



0J-A[B] cxp 



-B 



(10) 



with A — 3.8 and B = 0.94 [221 ■ Note that the validity of the WKB approximation is limited 
to S/S < —2, above which the zero point energy of a harmonic approximation in one of the 
energy minima reaches the barrier height. In practice the perfectly degenerate situation is 
hardly achieved, introducing some bias energy s which singles out a preferred ground state in 
which the beam is bent either to the left or right. We have thus an effective two state system 
with Hamiltonian 

H = -^a^-'-a, (11) 

and eigenenergies separated by -y/ (HA)^ + e"^, {ax and az are the Pauli spin matrices) . In the 
absence of any asymmetry s its eigenstates are coherent superpositions ^ \L) ±\R) of the az 
eigenstates \L) and \R), in which the nanotube is bent towards the left or right. The two states 
\L) or \R) are clearly macroscopically distinct [201 • Similar to experiments on flux qubits in 
SQUID rings |27[I28| . the existence of linear superpositions of these states may indirectly 
be verified by observing the avoided level crossing at e = 0. Such a macroscopic quantum 
coherence experiment with SWNT's requires that the small level splitting A can be detected 
against noise and damping in a mechanical resonance experiment and moreover, that the 
asymmetry e can be tuned through zero from any accidental nonzero value by external means. 
As shown in ^H], spectroscopy of the transverse vibrations of nanotubes is in principle possible 
by applying dc- plus ac-voltages on a charged nanotube. Moreover with a capacitive coupling 
the bias energy s may be changed via an appropriate electrostatic gate potential. Due to the 
still large mass involved, the tunnel splitting is rather small (around A = 27r • IMHz for 6/6 = 
—3), and thus coherent superpositions with an accessible value of A require nanotubes close 
to the buckling instability. Regarding the influence of damping effects, it is known ^29, that 
the dynamics of a two-level system subject to an ohmic dissipation mechanism is determined 
by the size of the parameter a = riq^/{2'KK). Here go = '^Urainl is the distance between the 
two minima and 77 the phenomenological damping parameter which appears in the equation 
of motion 

d 

m^sAi + riAi + -TTT-V^Ai) = 0. (12) 

Coherence in the two level system is present only for a < i at T = and ksT < hA/a for 
finite temperature and very small a. This requires that the quality factor of the SWNT in 
the uncompressed case (which is related to 77 by Q = 'meSi^o/v) obeys Q > 4\6\/{7r6^/'^). For 
the above value of 6, this leads to Q > 220 in the relevant regime \6\ ~ 6. This condition does 
not seem too stringent for SWNTs, note that a quality factor of Q = 500 has been reached 
for Si-based resonators in the GHz regime |2()j . 
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Conclusions. - We have discussed quantum effects in the mechanical properties of single 
wall carbon nanotubes, in particular zero point fluctuations in the transverse vibrations and 
the possibility to see the analog of MQC in nanobeams below the Euler buckling instability. 
While thermal vibrations of clamped SWNT's down to lengths L = 0.5/im have indeed been 
observed very recently : 16^ , it remains a considerable challenge to measure the tiny zero point 
vibrations of order 0.1 nm predicted for SWNTs of length L = 0.1/im near the buckling 
instability. With the sensitivity attained very recently with Si-based resonators [23 ! however, 
reaching this goal in the near future seems to be quite realistic. As regards the possibility to 
see the analogue of MQC near or below the Euler buckling instability, this requires to tune 
these systems rather closely below the instability point and doing spectroscopy with both dc- 
and ac-driving. Provided the methods used for multiwall nanotubes with lengths of several 
/xm |15j can be scaled down to clamped single wall nanotubes, quantum mechanics in its literal 
meaning would finally be of relevance in truly mechanical devices. 
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